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I. INTRODUCTION 

The equations of the Einsteinian cosmology are derived from General Relativity, with the 
use of the Cosmological Principle. This formulation is geometric, so that the descriptions of 
the cosmological effects are associated to the geometry of the spacetime, which is quantified 
by its curvature^. In this case, differential geometry plays a fundamental role together with a 
wide range of some analytical tools which arise from geometry. In this metric approach, the 
models of the universe are then classified in accordance with the features of the curvature 
tensor. 

In Newtonian cosmology, the equation that describes the evolution of the universe is 
obtained from the equation of motion and the energy equation for particles (galaxies) subject 
to gravitational forces and considering the Cosmological Principle as well^. This suggests 
that we can use the formalism developed by Lagrange and Hamilton, writing down the 
Lagrangian and Hamiltonian of the system in order to obtain the equation and integral of 
motion. 

The Newtonian cosmology, a model of the universe based on Newtonian mechanics, was 
well established in the third decade of last century^*^. In this context, the cosmological 
equations which describe the time evolution of the universe obtained by considering the 
Cosmological Principle, assuming that the pressure is zero and using Newtonian dynamics 
and gravitation. The results obtained in this case show that the universe could not be 
static, obeying an equation which algebraically is similar to Friedmann’s equation of Gen¬ 
eral Relativity, but obviously has different interpretation. This model based on Newtonian 
treatment provides a good description of the expanding universe in which pressure is zero^. 
When pressure is taken into account it is necessary to modify the basic equations of the 
pressureless case^ and adopt some assumptions^ in order to garantee the analogy between 
the Newtonian and Einsteinian approaches. 

It is worth call attention to the fact that to obtain the equivalence between Newtonian 
and Einsteinian cosmologies it is necessary to combine the equations of classical mechanics 
and the Cosmological Principle which nowdays is strongly supported by the smooth Cosmic 
Microwave Background Radiation (CMB) whose spectrum is isotropic and uniform to^*^ 
within (AT/T) < 10“^. Thus, the fact that this anisotropy agrees very well with that of 
the initially homogeneous and isotropic universe indicates that the Cosmological Principle 
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should be valid at least over length scales much bigger than hundreds of Megaparsecs. 

Quantum Newtonian cosmology is the theoretical framework within which it is possible 
to And the wave function which predicts the behavior of the Newtonian universe. This wave 
function for Newtonian cosmology is obtained in non-relativistic quantum mechanics, and 
thus, all we need is to solve the Schrodinger equation for the system under consideration!^^—. 

A Heun equation is a second-order linear differential equation with four regular singu¬ 
larities!^. The confluent form of the Heun equation is obtained when two singular points 
coalesce and the point at infinity becomes irregular. In this case, we have the biconfluent 
Heun equation!^'!^. The Heun functions have a wide range of applications in physics, from 
quantum mechanics to cosmology!^ and the number of articles using their solutions have 
been growing, even in a scenario in which the theory concerning these functions is not yet 
completed. 

This paper is organized as follows. In section 2 we obtain the Hamiltonian operator for 
the Newtonian universe. In section 3, we solve the Schrodinger equation in this background. 
We present the biconfluent Heun equation, and then we obtain the Heun polynomials in 
section 4. In section 5 we obtain the energy spectrum for a particle (galaxy) moving in a 
Newtonian universe. In section 6, we write the general and exact expression of the wave 
function. Finally, in section 7, we present our conclusions. 

II. HAMILTONIAN OPERATOR 

In our previous paper—, using some assumptions about the structure of the Newtonian 
universe, we obtained the Lagrangian, L, for the motion of a particle (galaxy) of mass m, 
which is written as 

L{R, R) = -mR? + + -AmR^ , (1) 

2 R 6 

where A is the cosmological constant, and the last term in Eq. ([T]) corresponds to a kind of 
cosmological potential energy. Using this Lagrangian, the equation for the scale parameter, 
R{t), can be written as 

R=~nGpR + hR, ( 2 ) 

where M = AnR^p/S is the mass of the sphere (total mass of the Newtonian universe). It 
is worth call attention to the fact that we are considering the case p = 0, that is, a universe 
in which the pressure is zero. 


3 



Quantum Newtonian cosmology and the biconfluent Heun functions 


The Lagrangian given by Eq. ([T]) does not depend explicitly on time, because the system 
is under the action of a uniform force held, so that the constant of motion, H^, called 
classical Hamiltonian of the system, can be defined as 

(3) 

From Eq. ([3]), we can obtain an equation for the scale factor, R{t), which is analogous to 
the Friedmann equation, whose form is 

= I + -k, (4) 

where C = SnGpR^/3 and k = —2E/m are constants, and E = Hd is the total energy. 
Thus, from Eq. ([3]) we can write the classical Hamiltonian as 

rr A\ ^ '9 GMUl 1 . , , , 

Hd{R, R) = -mR^ - - - -AmR^ . (5) 

/ R b 


Now, we can construct a phase space as {R,Pr), where R is the comoving coordinate 
with 

Pr = mk (6) 


being the linear momentum. Thus, substituting Eq. (I6]) into Eq. ([5]), the classical Hamilto¬ 
nian takes the form 


Hd{R,PR) = ^ 


GMm 

R 


-AmR^ . 
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Therefore, performing the canonical quantization 


(7) 


Pr Prop 



( 8 ) 


we obtain the Hamiltonian operator for a particle (galaxy) moving in the Newtonian universe, 
which is given by 


H 


(P _ GMm 
2m dP? R 
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III. SCHRODINGER EQUATION 


(9) 


In this section we will solve the Schrodinger equation for a particle (galaxy) in the New¬ 
tonian universe, given by 


H%Ij{R) = EP{R) , 


( 10 ) 
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where it is understood that the eigenfunction 'ijj{R) corresponds to the eigenvalue E, with 

being the general wave function of the time-dependent Schrodinger 


equation. 

Substituting Eq. 


into Eq. (Bill. we obtain 


2m dR? 


+ 


GMm 1 2 


ij{R) = Eij{R) 


( 11 ) 


Note that all eigenfunctions correspond to bound states of positive energy. It is more suitable 
to rewrite Eq. flTTll in terms of dimensionless quantities. To do this, we first introduce the 
dimensionless parameter 


where 


2E 

m ’ 



( 12 ) 

(13) 


Let us also use the dimensionless variable 


X = tR , 


(14) 


where the coefficient r is given by 


r = 


mQ\ 


h 


1/2 


■ 


Then, Eq. flTT]) turns into 


_ _5_ I V,(..) = 0 , 


where the parameter 5 is given by 


(5 = 


AGMm^ 

rh? 


( 15 ) 


( 16 ) 


(17) 


Now, let us analyse the asymptotic behaviour of ijj when |a;| —?■ oo. For any finite value 
of the total energy E, the quantity 7 and the term proportional to x~^ becomes negligible 
with respect to in the limit |x| —?■ 00 , so that in this limit, Eq. (1T6|) reduces to 

d^ 


dx'^ 


— x^ tjj{x) = 0 


The functions which satisfy this equation for large enough |a;| are 

'iIj{x) = . 


(18) 


(19) 
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The wave function must be bounded everywhere, including x = ±oo. Therefore, the 
physically acceptable solutions must contain only the minus sign in the exponent. The 
analysis of the asymptotic behaviour of the system under consideration suggests that we 
have to look for solutions of Eq. (1T6|) valid for all x and having the form 

i>{x) = X e~^^^‘^P{x) , (20) 


where P{x) is to be determined appropriatelly. Substituting Eq. 
conclude that P{x) must satisfy the equation 


(PP 


— 2x 


dP 


+ 


(-,-3 )-14 


dx'^ \x J dx 
whose analytical solutions will be obtained in the next Section. 


P = 0 , 


into Eq. flTOl) . we 


( 21 ) 


IV. BICONFLUENT HEUN EQUATION 


Equation fl?I|) is a particular case of the biconfluent Heun equation^, which in the canon¬ 
ical form is given by 


dx^ \ X ^ J dx 

+ I (7 — CK — 2) — -[5 -I- (1 -I- q;)/3] — I y = 0 , 

where y{x) = HeunB(Q;,/3, 7 , 5; x) are the biconfluent Heun functions. 
Let us assume that the solutions of Eq. fl2^ can be written as 




= 2_^CsX 
s=0 


and thus, we have 


d 




d^ 


00 

= ^ SCgX 

S = 1 

00 


s—1 


= 5^s(s-l)c,a; 


s-2 


s=2 


Substituting Eqs. 


and fl25|) into Eq. ([22]), we arrive at 


s(s — l)CsX^ ^ ^ — /3 — 2x J SCsX^ ^ 


s=2 


S = 1 


+ ^ (7 - « - 2 ) - -[(5 (1 + «) 4 ]- 


N 00 

I ^ CsX" = 0 , 

^ s=0 


( 22 ) 


( 23 ) 

( 24 ) 
( 26 ) 


( 26 ) 
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from which follows 


s(s - ^+ 


ajsCgX 


s-2 


l3sCsX^ ^ 


s=2 


S=1 


5 = 1 


2sCsx'' + y ^(7 - a - 2)csx'‘ - ^[<5 + (1 + a)/3]csX^ ^ = 0 , 


S = 1 


s=0 


5 = 0 


or, equivalently 


y^(S' + 2)[S + l)c5+2x‘^ + y~^ (1 + (y){S + 2)c5+23^'^ 

=0 s=-i 

OO CX) OD 

y^/3(S'+ l)cs+ix^ -'Y^2Scsx^ + y^(7 -a- 2)csx^ 


s=o 


s=o 

OO 


5=1 


5=0 


“ + oi)P]cs+ix^ = 0 , 


5=-l 


where S = s — 2. Collecting all terms of same order in x, we get 


(27) 


(28) 


(1 + a)ci - ^[(5 + (1 + a)/9]co \ x ^ 


+{2(2 + a)c 2 ~ ^ 2 + j3\ci + {'^ — a — 2)co} 


OO . 

y^{(S' + 2)(S' + 2 + a)c5+2 - ^ + (1 + 

5=1 


oi)(3\ + (S' + l)/3 > C 5+1 


+ (7 - a - 2 - 2 ^)cs}x^ = 0 . 


(29) 


Thus, Eq. (|29D gives the following recursion relations for the expansion coefficients 


(1 + a)ci — -[(5 + (1 + a)/3]co , 


2(2 + a)c 2 — 2 ^)P\ + /3f’Ci — (7 — a — 2 )co , 


(S' + 2)(S' + 2 + a)c5+2 — 2 + (S' + l)/3 j C 5+1 

-( 7 -a- 2 - 2 ^)c 5 , ^> 0 . 


(30) 
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Considering an appropriate choice of Cq, namely Cq = 1, relations given by Eq. fl30|) turn into 

Co = 1 , 

(1 + Q!)ci = -[5 + (1 + <y)l3] , 

2(1 + Q !)(2 + Q !) c 2 = + (1 + + /^ l ’ + (1 + *^)/^] 

-(1 + q ;)(7 - a - 2 ) , 

6(1 + tt)(2 + a)(3 + cr)c 3 = ^')P\ 2/3^ 

x{|+ (1 + ci;)/?] + /^l -[5 + (1 + «)/3] — (1 + a)(7 — a — 2)} 

—2(2 + Q !)(7 — a — 4)-[(5 + (1 + a)/?] . (31) 

Thus, taking into account a convenient change in notation, we have 

^0 = Co = 1 , 

^1 = (1 + cr)ci = 2 + ci)/3] , 

2^2 = 2(1 + cr )(2 + cr)c 2 , 

^3 = 6(1 + cr)(2 + cr)(3 + cr)c 3 , 


= s!(l + a)sCs, s > 0 , 


(32) 


where 


(1 + a)s 


r(s + 1 + tt) 
r(i + tt) 

j (1 + cr )(2 + cr) ... (s + «) 


s = l,2,3,... 
s = 0 . 


(33) 


Now, if we assume that a is not a negative integer, the biconfluent Heun functions can 
be written as^i 


HeunB(a, /3, 7 , 6; x) 


E 


As 

(1 + a)s s! 


(34) 


where 


As+2 — |(s + l)/3 + + *^)/^]| ^s+1 

— {s + l)(s + 1 + a )(7 — CK — 2 — 2s)As, s > 0 . 
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A. Heun polynomials 


From the recursion relation given by Eq. (I35|) . the function HeunB(a,/3, 7 , 5; x) becomes 
a polynomial of degree n if and only if the two following conditions are fulfilled^: 


(i) 7 — a — 2 = 2n, u = 0,1, 2, 

(ii) JU+i = 0 


(36) 


where A„_|_i is a polynomial of degree n + 1 in 5. Note that there are at most n + 1 suitable 
values of 5, labeled 

5 ;, 0<p<n. (37) 

The polynomial which has n + 1 real roots when 1 + a > 0 and /5 G M, is the 

determinant of dimension n + 1 given by 


5' 1 0 0 ... 

2(1 + a)n 5' — (3 1 0 

0 4(2 + a)(n-l) 5'- 2(3 1 0 


0 


0 


72 5' -3(3 1 


0 

0 

0 

0 


0 


: 7s-i 1 

0 0 js 


(38) 


where 


S' = --[(5 + (1 + a) 


(39) 

(40) 

(41) 


5( = <5'-(s + l)/3, 

7 s = 2(s + l)(s + 1 + a){n - s) . 

Now, let us return to the wave functions '0(x). Using Eq. fl20|) . we see that the physically 
acceptable solutions of Eq. (1^ are given by 

(42) 




where the functions (3] x) are polynomials of degree n satisfying the biconfluent Heun 

equation given by Eq. fl2^ . with 7 = 0 ; + 2(n + 1), A„+i = 0, and 0 < p < n, namely: 






dx"^ 


+ 


1 T CK 


X 


- (3-2x 


dP, 




dx 


+ 1 2n — -[(5” + (1 + o)(3] — > Pn,ij. — 0 


X 


(43) 
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The polinomials /9; x) = HeunB(a, /?, a + 2(n + 1), 5”; x) are called Heun polyno¬ 

mials of the biconfluent case. It is clear from the foregoing discussion that they are uniquely 
defined, except for an arbitrary multiplicative constant. Its highest order terms are given 
by 




- 6 ’ 


X 


n—l 


+ 


-[(5Q5' — /3) — 2(1 -|- a)n] 


X 


n—2 


(1 + Of) 


(l-Fa)i" 1! ' (1a): 

^ [P((5' - /3) {6' - 2/3) - 4(5'(2 + a)(n - 1) 


2 ! 


X 


n—3 


-2{6' — 2/3)(l -|- a)n]——h Q{x) , 


(44) 


where Q{x) is a polynomial with degree < u — 4. 

In our case, the parameters a and /3 take the values 1 and 0, respectively. Therefore, the 
explicit form of the Heun polynomials P„_^(l, 0; x), for n = 0,1, 2, are as follows: 


• n = 0, 7 = 3, (fg = 0 


Po,o(l,0;x) = 1 ; 


(45) 


• u = 1, 7 = 5, ((f^)^ — 16 = 0 ^ (fg = +4, = —4 

(46) 

(47) 

. n = 2, 7 = 7, {6lf - 8061 = 0^6^ = +4^5,6f = 0,<5| = -4v^ 


Pi,o(l,0;x) = x-]~ 


Pi,i(l,0;x) =x- 


2 2 


X — 1 , 

X -h 1 ; 


. , 9 1 (5n 1 

P 2 0 ( 1 , 0; x) = x^-^x -I- - 

Auv > j y 226 


(<5. 


2\2 

0 / 


P 2 ,i(l, 0 ;x) = 


X 


2 UJ 1 
2 2'"'+6 




. ^ 9 1 <^2 1 

^ 2 , 2 ( 1 , 0 ;x) = x^ - 2 y^+ 6 


( 5 , 


2\2 


8 


= x^ — \/5x -I- 1 , 


2 - 3 ’ 


- = x^ -I- \/5x -I- 1 . 


(48) 

(49) 

(50) 


It is worth noting that for all values of n, namely, n = 0,1, 2,..., correspond a function 
HeunB(a, /3, 7 , 5; x) which is a polynomial of degree n in x. 
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FIG. 1. The effective potential energy of the Newtonian universe Ueff{R) = — + ||A|mi?^ 

and the first five energy levels are shown. 

V. ENERGY LEVELS 

In order to obtain the energy spectrum, let us use Eqs. (|T^ and (j36l) . and take into 
account that a = 1. Thus, the energy spectrum for a particle (galaxy) in the Newtonian 
universe is given by 

I 

u+0 hfl = (n+0 h (-^) , n = 0,1,2,.... (51) 

We see from Eq. fl5T]) that this quantum mechanical energy spectrum consist of an infinite 
sequence of discrete levels (see Fig. [T]), which are equally spaced and are similar to those 
obtained for the three-dimensional isotropic harmonic oscillator. Note that the eigenvalues 
given by Eq. flHTll are degenerate. 

If we consider a scenario in which the cosmological constant is negative, that is, A = — |A|, 
we can rewritten the energy spectrum as 

0 h , n = 0,l,2,..., (52) 

where (|A|/3)^/^ plays the role of the frequency if we compare Eq. (15^ with the corresponding 
to the energy spectrum of the three-dimensional isotropic harmonic oscillator. 
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FIG. 2. The first six wave functions of the Newtonian universe. 

VI. THE WAVE FUNCTION FOR NEWTONIAN COSMOLOGY 

Using Eq. fH2|) . we see that to each discrete value given by Eq. fl5T]) . there corresponds 
n + 1 physically acceptable eigenfunctions (see Fig. |2]), given by 

= Nn tR PnA(^A-,rR) 

= Nn tR HeunB(a, /3, a + 2(n + 1), tR) , (53) 

where the parameters a, A li 5 are given by the following values and expressions: 

a = 1 ; (54) 

/9 = 0 ; (55) 

(- 1 )“'*. ( 57 ) 

The quantity A„ in Eq. (I53|) is a constant which (apart from a arbitrary phase factor) 

can be determined by requiring that the wave function given by Eq. fl53l) obeys the following 
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relation 


\^l:nAR)\"dR = 


\Nr. 


12 poo 


X 


T 


“ PlAa, l};x)dx=l 


( 58 ) 


In order to evaluate this integral, we consider the power series expansion of Pn,^{a, (3;x) 
in X given by Eq. (I34|) . as well as the second power series expansion 




Using Eqs. 


where 


Ai x^ 

^ (l + a)z l\ 

and fl59|) . we obtain the following result 

/r\ 1/2 

Nn = (y) , 

poo 

1=1 x‘^e~''^Pn,fx{a,l3]x)Pm,u{a,^-,x)dx 

Jo 

"" "" A.A, 1 /•°° 


( 59 ) 


( 60 ) 


EE 


^ i^Q (1 + «)fc(l + kill Jq 


X 


2+k+l^-x 


(61) 


For n = m, fi = u and k = I, we have 


/ = 


X e ^ P {a, ^■,x)dx 


E 


Al 


[{l + a)kk\Y 


x2+2k^-x^dx 


k=0 


Since 


the integral in Eq. 


'0 


u \u 


is simply 

^ = E 


Al 


Co 1(1 


+ k 


( 62 ) 


( 63 ) 


( 64 ) 


From Eqs. flbOp and fl64p we see that apart from an arbitrary complex multiplicative factor 
of modulus one the normalisation constant i\Q is given by 

- 1/2 


N„, = 


E 

.k=0 


Al 


[(1 + a)fcA;!]2 2r 


r hr + fc 


(65) 


Using Eq. (165]), we can write the normalized solution of the Schrodinger equation in a 
Newtonian universe as 


(-R) 


Al 


2 I(l + a)fc/d]2 2r^ 


k 


- 1/2 


X tR e HeunB(Q!, /3, a + 2(?7, + 1), tR) . 


( 66 ) 
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VII. CONCLUSIONS 

We have presented exact and general solution of the Schrodinger equation for a particle 
(galaxy) moving in a Newtonian universe in the presence of a cosmological constant term. 
The complete set of solutions, given in terms of the biconfluent Heun functions, satifles 
the appropriate boundary conditions -0(0) = 0 and = 0, and is valid over the range 

0 < a: < cxo. 

In order to get a well-behaved solution, we imposed the polynomial condition on these 
analytic solutions and thus we obtained the Heun polynomials for the biconfluent case. 

An exact expression for the energy spectrum was also obtained. These energy levels are 
similar to the ones corresponding to a three-dimensional isotropic harmonic oscillator, with 
an effective frequency given in terms of the cosmological constant. 

As a conclusion we emphasize that the wave functions as well as the energy spectrum 
codifies the presence of a global cosmological force which affects the particles (galaxies) of 
this Newtonian universe. 
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